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ABSTRACT 

The  reversal  mechanism  for  the  magnetization  in  an  Fe  film  with  an  array  of  elliptical  holes  is  investigated  using  the  diffracted 
magneto-optic  Kerr  effect  (D-MOKE)  technique.  D-MOKE  results  are  obtained  as  a 

function  of  temperature  and  the  angle  between  the  applied  magnetic  field  and  an  ellipse  axis.  The  transverse  and  longitudinal  magnetization 
components  and  minor  magnetization  loops  are  also  explored  in  order  to  understand  the  reversal  process.  The  experimental  results  are 
interpreted  using  micromagnetic  simulations.  The  simulations  account  for  the  strong  angular  dependence  of  the  hysteresis  loops  and  provide 
a  detailed  picture  of  how  the  local  magnetization  evolves  during  reversal.  The  actual  reversal  process  occurs  neither  by  coherent 
rotation  of  domains  nor  by  clear  domain-wall  motion:  domain  smearing  appears  to  be  a  more  suitable  description  of  the  phenomenon. 
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The  reversal  mechanism  for  the  magnetization  in  an  Fe  him  with  an  array  of  elliptical  holes  is  investigated 
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and  longitudinal  magnetization  components  and  minor  magnetization  loops  are  also  explored  in  order  to 
understand  the  reversal  process.  The  experimental  results  are  interpreted  using  micromagnetic  simulations.  The 
simulations  account  for  the  strong  angular  dependence  of  the  hysteresis  loops  and  provide  a  detailed  picture  of 
how  the  local  magnetization  evolves  during  reversal.  The  actual  reversal  process  occurs  neither  by  coherent 
rotation  of  domains  nor  by  clear  domain-wall  motion:  domain  smearing  appears  to  be  a  more  suitable  descrip¬ 
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INTRODUCTION 

Arrays  of  micron-  and  submicron-sized  holes  in  a  con¬ 
tinuous  ferromagnetic  film  have  been  proposed  as  a  competi¬ 
tor  for  high-density  storage  media,  with  characteristics  of 
high  stability  while  avoiding  the  superparamagnetic  limit.1 
The  introduction  of  arrays  of  holes  in  magnetic  thin  films 
also  provides  a  means  of  engineering  their  magnetic  proper¬ 
ties  in  a  controllable  way.1'2  The  switching  mechanism  dur¬ 
ing  magnetization  reversal  is  an  important  issue,  which  is  not 
yet  well  understood  in  these  systems. U6  In  a  recent  article 
we  presented  the  results  of  a  diffracted  magneto-optic  Kerr 
effect  (D-MOKE)  investigation  of  an  array  of  elliptical  holes 
in  an  Fe  film.7  There  we  showed  how  the  D-MOKE  signal  is 
related,  through  the  magnetic  form  factor,  to  the  domain 
structure  that  exists  during  reversal.  We  also  showed  that 
when  the  field  is  applied  along  the  long  axis  of  the  ellipses 
only  a  small  fraction  of  the  sample  develops  domains  while 
extensive  “stripe”  domains,  that  connect  next-nearest  el¬ 
lipses,  form  when  the  field  is  applied  along  the  ellipse’s  short 
axis.  In  that  investigation  no  attempt  was  made  to  extract 
from  the  D-MOKE  results  how  the  domains  formed  and/or 
transformed  during  reversal. 

The  sample  investigated  here  is  the  same  one  investigated 
in  Ref.  7,  viz.,  an  Fe  film  with  nanometric  elliptical  holes 


forming  a  square  1  X  l-/zm2  lattice.  Here  we  report  a  com¬ 
prehensive  study  of  the  D-MOKE  results  as  a  function  of  the 
angle  between  the  applied  field  and  the  ellipse  axis,  the  tem¬ 
perature,  the  transverse  and  longitudinal  magnetization,  and 
minor  magnetization  loops  (i.e.,  loops  in  which  the  maxi¬ 
mum  and  minimum  fields  are  not  equal).  The  results  are 
interpreted  using  the  form  factor  approach  described  in  Ref. 
7.  The  analysis  of  D-MOKE  loops  has  been  improved  by 
extracting  field-dependent  form  factors  from  micromagnetic 
simulations.  In  spite  of  the  difficulties  associated  with  micro- 
magnetic  simulations  of  negative  arrays  (viz.,  dealing  with 
the  boundary  conditions),  such  calculations  enable  full  hys¬ 
teresis  loops  from  any  interference  order  to  be  calculated, 
and  simultaneously  yield  a  detailed  picture  of  how  the  local 
magnetization  evolves  during  reversal. 

EXPERIMENTAL  DETAILS  AND  PREVIOUS  RESULTS 

Details  of  the  sample  fabrication  were  previously  de¬ 
scribed  in  Ref.  7.  Briefly,  the  sample  is  an  Fe  film,  60  nm 
thick,  with  elliptical  holes  ~-200  nm  wide  by  «*800  nm  long, 
forming  a  square  1  X  1  -/zm2  lattice.  A  transmission  scanning 
electron  microscopy  image  of  the  sample  is  shown  in  Fig. 
1(a). 

MOKE  experiments  were  performed  in  the  “transverse 
Kerr”  geometry.  Incident,  p-polarized  light,  on  the  sample 
was  detected  in  the  scattering  plane  with  no  analyzer.  This 
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FIG.  1.  (a)  Scanning  transmission  electron  microscopy  (STEM) 
image  of  our  array  of  elliptical  hoes  in  a  60-nm-thick  Fe  film,  (b) 
Lorentz  STEM  image  of  the  same  sample;  the  arrows  are  the  mag¬ 
netization  directions  inferred  from  D-MOKE  analysis. 

configuration  is  sensitive  only  to  the  component  of  M  per¬ 
pendicular  to  the  scattering  plane.  By  applying  a  field  (//) 
perpendicular  to  the  scattering  plane  we  obtain  the  magneti¬ 
zation  parallel  to  H  (longitudinal  magnetization  component), 
with  H  in  the  film  plane  and  in  the  scattering  plane  we  obtain 
the  component  of  M  perpendicular  to  H  (transverse  magne¬ 
tization  component). 


The  basic  formalism  developed  to  interpret  D-MOKE  re¬ 
sults  was  presented  in  Refs.  7  and  8.  In  Ref.  7  the  domain 
structure  at  remanence  was  inferred  from  simple  energy  con¬ 
siderations.  In  Ref.  8  and  here  we  obtain  form  factors  as  a 
function  of  field  from  micromagnetic  simulations.  Micro- 
magnetic  calculations  were  performed  using  the  web  version 
of  the  NIST  code.9  The  difficulties  of  using  this  approach  for 
negative  arrays  will  be  discussed. 

Figure  1(b)  is  the  Lorentz  microscopy  image  showing  the 
domain  structure,  deduced  from  the  D-MOKE,7  that  forms 
when  the  field  is  applied  along  the  short  ellipse  axis.  The 
form  factors  of  these  domains  provide  a  qualitative  explana¬ 
tion  of  the  salient  features  of  the  longitudinal  D-MOKE 
loops  shown  in  Figs.  2(d)-2(f).  Because  the  domains  that 
form  when  the  field  is  applied  along  the  long  axis  of  the 
ellipses  are  quite  small,  their  effect  on  the  D-MOKE  loops 
[Figs.  2(a)-2(c)]  is  also  small.  For  this  reason,  in  this  paper, 
we  will  concentrate  on  the  magnetization  reversal  for  fields 
that  are  roughly  along  the  ellipse’s  short  axis. 

RESULTS 

Angular  dependence  of  transverse  and  longitudinal 
magnetization 

Figure  3  shows  the  zeroth-  and  first-order  longitudinal  and 
first-order  transverse  loops  measured  as  a  function  of  the 
angle  ( (/>)  between  the  applied  field  and  the  ellipse’s  short 
axis.  (The  zeroth-order  transverse  loops  are  noisy  but  mimic 
the  first-order  ones).  For  (f>  =  0°  ±  1.5°  we  find  that  the  loops 


FIG.  2.  Longitudinal  D-MOKE  loops  of  vari¬ 
ous  orders,  from  the  patterned  area  of  the  film. 
The  applied  field  is  along  the  long  and  short  axes 
of  the  ellipses. 
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FIG.  3.  Experimental  D-MOKE  loops  vs 
angle  {(/>)  between  field  and  short  ellipse  axis: 
zeroth-  and  first-order  longitudinal  and  first-order 
transverse. 


have  quite  different  shapes  than  for  <j>>2° .  As  will  be  shown 
in  a  quantitative  fashion  in  the  next  section,  these  differences 
can  be  traced  to  the  manner  in  which  “domains”  are  initially 
nucleated.  Qualitatively,  as  discussed  in  Ref.  7,  the  increase 
in  the  first-order  longitudinal  D-MOKE  signal  at  low  values 
of  H  with  respect  to  its  value  when  the  sample  is  saturated 
(high  values  of  H ),  is  due  to  the  formation  of  domains.  The 
different  shapes  of  the  loops  in  Fig.  3  indicate  that  when  the 
field  is  along  the  hard  axis  (<fi=  ±  1.5°),  the  onset  of  mag¬ 
netization  reversal  is  retarded,  and  considerably  more  do¬ 
main  formation  occurs  prior  to  reversal.  As  the  angle  <f>  is 
varied  (either  increased  or  decreased),  domains  are  present 
both  before  and  after  reversal,  and  the  value  of  H  at  which 
the  reversal  occurs  becomes  lower  as  cj)  increases. 

Minor  loops 

Loops  in  the  left  column  of  Fig.  4  show  experimental 
first-order  longitudinal  D-MOKE  loops  as  a  function  of  the 
maximum  positive  field  applied  during  the  hysteresis  cycle. 
While  the  zeroth-order  loops  (not  shown)  show  only  minor 
changes  as  the  maximum  field  of  the  loop  is  changed,  the 
first-order  loops  show  substantial  differences  close  to  the  re¬ 
versal  field.  In  particular,  the  negative  overshoot  on  the  de¬ 
creasing  side  of  the  loop,  which  is  almost  absent  in  the  full 
loop,  becomes  a  prominent  feature  for  the  asymmetric  loops. 
These  changes  will  again  be  traced  to  the  way  in  which  do¬ 
mains  are  nucleated  and  annihilated. 


Temperature  dependence 

Figure  5  shows  zeroth-  and  first-order  longitudinal 
D-MOKE  loops  as  a  function  of  temperature  ( H  is  applied  at 
a  small  angle  from  the  short  axis  of  the  ellipses).  We  observe 
no  substantial  changes  in  the  hysteresis  loops  as  the  tempera¬ 
ture  is  lowered.  As  discussed  in  Ref.  7,  the  magnetization 
reversal  occurs  via  the  nucleation  of  demagnetizing-field- 
induced  90°  domains,  which  form  to  reduce  the  dipolar  mag¬ 
netostatic  energy  associated  with  the  surface  charges  that  ap¬ 
pear  at  the  hole  edges.1011  The  minor  modifications  of  the 
D-MOKE  loops  vs  temperature  provides  evidence  that  ther¬ 
mal  activation  does  not  play  a  dominant  role  in  the  domain 
formation  during  reversal. 

DISCUSSION 

Micromagnetic  simulations  and  Theory  of  D-MOKE 

We  remind  the  reader  that,  as  described  in  Ref.  7,  the 
D-MOKE  signals  are  proportional  to  the  magnetic  form  fac¬ 
tors  if)  defined  by 

f=  j  myexp[inG’r]dS,  (1) 

where  my  is  the  magnetization  per  unit  area  along  the  direc¬ 
tion  perpendicular  to  the  scattering  plane,  n  is  the  interfer¬ 
ence  order,  G  is  the  reciprocal-lattice  wave  vector,  and  the 
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FIG.  4.  Experimental  (left)  and  simulated 
(right)  first-order  longitudinal  loops  as  a  function 
of  maximum  applied  positive  field. 


integration  is  carried  out  over  a  unit  cell  of  the  array.  Note 
that  for  n  =  0  (i.e.,  the  reflected  spot)  Eq.  (1)  is  just  the  inte¬ 
gral  of  My,  and  hence  yields  the  same  information  as  con¬ 
ventional  MOKE.  In  Ref.  7  it  was  assumed  that  at  some 
instant  during  each  hysteresis  loop  domain  formation  was 
complete  in  the  sense  that  it  was  the  same  in  all  unit  cells. 
With  this  assumption  it  was  shown  that  the  domains  dis¬ 
played  in  Fig.  1(b)  were  consistent  with  the  maximum  in¬ 
creases  or  decreases  observed  in  the  loops.  No  attempt  was 


made  in  that  investigation  to  explain  how  domains  appear 
and  evolve  during  the  switching.  Here  we  wish  to  examine  if 
the  path  taken  during  domain  formation  can  be  extracted 
from  the  data  presented  in  Figs.  2-4. 

The  field-dependent  form  factors  can,  in  principle  be  ex¬ 
tracted  from  micromagnetic  simulations.8  At  each  field  the 
distribution  of  spins  is  extracted  from  the  simulation  and, 
using  Eq.  (1),  the  form  factors  obtained.  Although  this  ap¬ 
proach  works  remarkably  well  for  a  system  of  circular 
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FIG.  5.  Temperature  dependence  of  zeroth- 
and  first-order  longitudinal  loops. 
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FIG.  6.  Calculated  loops:  zeroth-  and  first- 
order  longitudinal,  and  first-order  transverse  for 
4>  equal  to  1°  and  10°. 


disks,12  it  encounters  some  difficulties  when  applied  to  an 
array  of  holes  as  in  the  present  case.  The  origin  of  the  diffi¬ 
culty  lies  in  establishing  appropriate  boundary  conditions 
(BCs).  From  a  calculation  standpoint  it  is  desirable  to  per¬ 
form  the  simulation  on  a  single  unit  cell  of  the  negative 
array.  This  approach  requires  the  introduction  of  periodic 
BCs  necessary  to  deal  with  an  extended  structure.  Introduc¬ 
ing  the  correct  BCs  for  this  approach  to  be  valid  is,  however, 
far  from  trivial.  Conceptually  it  would  be  easy  to  introduce 
the  exchange  coupling  at  the  boundaries  of  the  unit  cell  by 
coupling  each  edge  to  its  opposite  side.  The  dipolar  terms, 
however,  are  not  amenable  to  such  a  simple  solution  since 
they  depend  on  the  distribution  of  charges  throughout  the 
neighboring  cells.  In  the  absence  of  a  suitable  way  of  intro¬ 
ducing  the  BCs,  we  have  resorted  to  the  brute  force  approach 
in  which  a  large  system,  comprised  of  many  unit  cells,  is 
simulated.  Here  again  some  compromises  must  be  reached. 
Edge  effects,  that  enhance  magnetization  rotation  when  M  is 
perpendicular  to  an  edge  (i.e.,  favoring  domain  nucleation) 
and  hinder  its  rotation  when  M  is  parallel  to  an  edge  (i.e., 
inhibiting  domain  nucleation  or  rotation),  are  reduced  by 
simulating  a  structure  with  rough  edges.  The  roughening  is 
accomplished  by  alternately  removing  a  computational  cube 
from  around  the  edges.  Although  this  reduces  the  effect  of 


the  outer  edges  on  the  central  cells,  the  simulations  invari¬ 
ably  show  a  non  uniform  formation  of  domains  in  different 
unit  cells.  Ideally  the  size  of  the  system  could  be  increased 
until  a  uniform  domain  structure  is  found  for  the  central 
cells.  In  our  case,  even  for  a  5  X5  unit  cells  system  we  were 
not  in  that  regime.  Larger  systems  required  unacceptable  cal¬ 
culation  times.  Our  simulated  structure  consisted  of  5  X  5 
unit  cells  each  one  consisting  of  21X21  calculation  cells  (for 
a  total  of  1 1  025  computational  cubic  cells  each  one  with  a 
side  of  48  nm)  and  the  outer  edges  were  roughened.  The 
material  parameters  of  Fe,  with  the  anisotropy  set  to  zero, 
were  used  for  the  simulations.  The  form  factors  were  calcu¬ 
lated  for  the  central  unit  cell  of  the  5X5  simulation. 

Comparison  between  measurements  and  theory 

Figure  6  shows  the  zeroth-  and  first-order  longitudinal 
loops  and  first-order  transverse  loops,  calculated  from  the 
form  factors.  The  applied  field  was  at  1°  and  10°  from  the 
ellipse’s  short  axis.  The  loops  in  Fig.  6  qualitatively  repro¬ 
duce  the  experimental  hysteresis  loops  in  Fig.  3.  Both  experi¬ 
mental  and  simulated  loops  show  that  the  zero-order  loops 
are  not  greatly  affected  by  the  changes  in  field  direction.  For 
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FIG.  7.  Micromagnetic  domains  vs  H  (2000,  500,  and  — 100  Oe) 
and  angle  <f>  ( 1°  and  10°).  In  panels  (c)  and  (d)  a  gray-scale  propor¬ 
tional  to  the  magnetization  component  perpendicular  to  the  short 
axis  of  the  ellipses  is  used  to  highlight  the  domain  structure.  The 
dotted  lines  are  a  guide  to  the  eye. 

the  first-order  longitudinal  and  transverse  loops  the  shape 
changes  can  be  ascribed  to  variations  in  the  domain  structure 
during  reversal. 

In  Fig.  7  we  show  the  micromagnetic  spin  configurations 
at  various  fields  for  the  two  directions  of  the  applied  field. 
Figure  7(d)  shows  the  previously  reported7  stripelike  do¬ 
mains,  with  the  magnetization  vector  at  '“90'J  with  respect  to 
the  initial  saturation  direction.  These  domains  bridge  next- 
nearest-neighboring  holes  along  the  lattice  diagonal,  when  77 
is  applied  at  10°.  For  77  applied  almost  along  the  short  ellipse 
axis.  Fig.  7(c)  shows  that  domain  nucleation  is  hindered  and, 
when  it  does  occur,  is  different  from  that  for  the  10°  angle. 
Indeed,  the  comparison  of  Figs.  7(c)  and  7(d)  shows  that  for 
77  =  500  Oe  the  90°  domains  are  less  extended  in  the  case  of 
77  applied  almost  along  the  short  ellipse  axis,  and  they  form 
short  blades  that  do  not  overlap  with  those  from  next-nearest 
neighboring  holes.  Figures  7(e)  and  7(f)  show  that  when  77  is 
reversed  (77=  — 100  Oe)  an  intricate  domain  structure  is 
generated.  A  detailed  comparison  of  these  two  figures  shows 
that  the  average  magnetization  component  parallel  to  77  has 
changed  sign  (i.e.,  the  jump  in  the  hysteresis  loops  has  al¬ 
ready  occurred)  only  when  77  is  applied  at  10°.  This  is  in 
agreement  with  the  retarded  magnetization  reversal  observed 
experimentally  when  77  is  close  to  the  hard  axis.  Also  evident 
in  Fig.  7  is  that  reversal  does  not  occur  via  an  obvious  do¬ 
main  wall  motion  nor  by  coherent  rotation.  It  is  perhaps 
more  suitable  to  describe  the  reversal  process  as  domain 
smearing. 

The  experimental  full  loop  in  Fig.  4  is  indicative  that  if 
the  sample  approaches  reversal  after  being  saturated,  most  of 
the  ‘domain’  formation  occurs  before  reversal  (viz.,  the  peaks 
in  the  first-order  loops  are  larger  before  than  after  reversal). 
As  the  maximum  field  is  reduced,  the  intensity  of  the  two 
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FIG.  8.  Domains  in  minor  loops  at  the  same  field  (  —  300  Oe)  but 
for  different  max  fields:  panel  (a)  maximum  field  2500  Oe,  panel 
(b)  maximum  field  300  Oe. 

peaks,  after  and  before  reversal,  become  comparable  (middle 
loop).  Loops  in  the  right  column  of  Fig.  4  show  the  first- 
order  longitudinal  loops  calculated  using  the  form  factor  fx 
extracted  from  micromagnetic  simulations  as  a  function  of 
the  maximum  positive  field  applied  during  the  hysteresis 
cycle.  The  qualitative  agreement  with  the  measured 
D-MOKE  loops  is  remarkable.  In  Fig.  8  we  show  two  spin 
distributions  at  the  same  field  77  =—  0.3  kOe  but  reached 
from  different  maximum  fields:  2.5  [Fig.  8(a)]  and  0.3  kOe 
[Fig.  8(b)],  corresponding  to  the  loops  shown  in  Figs.  4(a) 
and  4(c)  respectively.  Also  indicated  in  Fig.  8  are  the  first- 
order  form  factors  for  the  two  states,  normalized  to  the  value 
of  fx  at  saturation  (we  actually  used  the  value  of  /]  at  77 
=  2.5  kOe.)  The  increase  of  fx  for  the  state  reached  from  a 
lower  maximum  field  [Fig.  8(b)]  reflects  the  larger  deviation 
from  the  saturated  state  and  provides  the  explanation  for  the 
loops  in  Fig.  4. 

CONCLUSIONS 

We  have  used  the  D-MOKE  technique  to  investigate  the 
reversal  mechanism  in  an  array  of  elliptical  holes  in  an  Fe 
film.  Experimentally,  based  on  the  changes  in  the  various 
diffracted  hysteresis  loops,  we  find  that  the  reversal  process 
depends  quite  strongly  on  the  angle  between  the  magnetic 
field  and  the  ellipse  axis.  In  spite  of  boundary  condition  dif¬ 
ficulties,  we  have  used  micromagnetic  simulations  to  inter- 
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pret  the  experimental  results.  The  micromagnetic  boundary 
condition  problems  were  partially  circumvented  by  introduc¬ 
ing  rough  edges  on  the  simulated  sample  and  by  using  the 
largest  sample  consistent  with  acceptable  computing  times. 
The  simulations  account  for  the  strong  angular  dependence 
of  the  hysteresis  loops.  They  also  show  the  angular  depen¬ 
dence  of  the  domain  structure  generated  during  reversal.  The 
actual  reversal  process  occurs  neither  by  a  coherent  rotation 
of  domains  nor  by  clear  domain  wall  motion:  domain  smear¬ 
ing  appears  to  be  a  more  suitable  description  of  the  phenom¬ 
enon. 
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